Multi-mode dynamics of a coupled ultracold atomic-molecular system 
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We analyze the coherent multi-mode dynamics of a system 
of coupled atomic and molecular Bose gases. Starting from an 
atomic Bose-Einstein condensate with a small thermal com- 
ponent, we observe a complete depletion of the atomic and 
molecular condensate modes on a short time scale due to sig- 
nificant population of excited states. Giant coherent oscilla- 
tions between the two condensates for typical parameters are 
almost completely suppressed. Our results cast serious doubts 
on the common use of the 2-mode model for description of 
coupled ultracold atomic-molecular systems and should be 
considered when planning future experiments with ultracold 
molecules. 



After experimental achievement of quantum degener- 
acy in atomic gases of bosons |ij and fermions |2| , sev- 
eral leading groups started the quest for creation of an 
analogous state in a gas of ultracold molecules H . Such 
systems present a potentially much more challenging ex- 
perimental task as there are several new degrees of free- 
dom to be controlled. A major step towards this goal 
has been made by the Austin group of Heinzen Q] , who 
photoassociated atoms in an atomic condensate into a 
selected internal molecular state producing for the first 
time molecules in the nK regime. Cold molecules are 
believed to be created as a transient state in experi- 
ments with Feshbach resonances as well || . In this novel 
area there are only a few theoretical works analyzing the 
process of making ultracold molecules in atomic conden- 
sates and predicting several interesting phenomena in the 
mixed atomic-molecular systems. Javanainen et al. in a 
series of papers || analyzed the efficiency of photoassoci- 
ation of an atomic condensate into its molecular counter- 
part using various theoretical schemes. Others predicted 
that the ground state of the hybrid system would have 
a soliton-like nature with liquid-like properties per- 
forming Josephson-like oscillations in response to a sud- 
den variation of a magnetic field Most frequently, 
coherent oscillations between the atomic and molecular 
condensates are envisaged P-|lTf . All the approaches em- 
ployed so far have used a 2-mode approximation (only 
the condensates involved) to describe the dynamics of a 
coupled ultracold atomic-molecular system. 

In this Letter we demonstrate that the 2-mode ap- 
proach is inadequate as a description of current exper- 
iments on stimulated production of cold molecules in 
atomic Bose-Einstein condensates. For typical param- 



eters the predicted oscillations between the two conden- 
sates are strongly damped due to significant population 
of excited atomic and molecular modes leading to a com- 
plete depletion of the initial condensate on a short time 
scale. The method we use, quantum-optical in spirit, can 
be regarded as a generalization of the classic Bogolubov 
approximation p"2| |. 

The second-quantized Hamiltionian for the hybrid 
atomic-molecular system confined in a box with peri- 
odic boundary conditions may be written in the following 
form: 
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where $ and 4" are atomic and molecular field opera- 
tors, respectively, V — L 3 is a volume of the system 
(L being a size of the box), f2 parameterizes a coupling 
between the two fields and g = characterizes the 

atom-atom interactions in the low-energy, s-wave approx- 
imation (a being the scattering length and m - the mass 
of the atom; in fact, a is slightly changed in external, 
e.g. optical, fields - in present experiments this correc- 
tion is small though |I(J ) . A similar Hamiltonian is used 
to describe a process of second-harmonic generation in 
nonlinear optics |l3|] . The atom-molecule and molecule- 
molecule collisions are not included as their parameters 
in the low-energy regime are unknown. However, as can 
be seen from the construction of the method, their in- 
corporation could be easily accomplished. Note that in 
this model a molecule is created when the positions of 
two atoms exactly coincide (for an example of a finite- 
range coupling approach see The fields $ and W 
are expanded in natural modes of the system - the plane 
waves: 
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where ak and 6k are bosonic annihilation operators for 
atoms and molecules, respectively, and k = with 
rii = 0,±1,±2,... (i — x,y,z). With this substitution 
the Hamiltonian assumes its final form: 
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where £ = 2m (x") 2 - After elimination of a fast 
time dependence with the substitution {a, 6} k = 
exp(— i£n 2 i){a, /3} k , the Heisenberg equations of motion 
for the operators a k and /3 k acquire the following form: 

a k = -2ffi^exp(-^|n- n'| 2 t)/3 k+k /a k , - (4) 
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The first hint about limitations of the 2-mode model 
comes from the following argument. As initially only the 
atomic condensate (the k = mode) is populated, one 
might naively suspect that the coupling would primarily 
lead to an interconversion between atoms and molecules 
H@liH0]l- ln sucn a case only the atomic and molec- 
ular condensate (k = 0) modes would be macroscopi- 
cally populated and therefore we replace the correspond- 
ing operators («o and Po) by c- numbers (classical com- 
plex fields) and set all the others to zero (a k ^o = and 
ftk^o = 0). Then, in the absence of atomic collisions, 
time evolution of the amplitudes can be calculated ana- 
lytically. Assuming ao = rexp(i^) and (3q = pexp(i#) 
with 4> — 4>q = const and 9 = 9q = const, the solution is: 
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where TV, the total number of atoms, is a conserved quan- 
tity (N = \a (t)\ 2 + 2\f3 (t)\ 2 ) - see Fig.g. Note that the 
inclusion of interactions makes the 2-mode problem ana- 
lytically insoluble (in particular the condition of constant 
phases 9 and <p cannot be fulfilled any more). Numeri- 
cal solutions are of oscillatory character [P-|i"lt - different 
from (^J) . They depict an interconversion between atomic 
and molecular condensate modes. 

In the next step we calculate quantum corrections to 
such a 2-mode model treating the condensates (ao and 
/?o) as the only sources of particles. This amounts to 
leaving in Eqs.([|) the terms with at least one k = (in 
the form of (||)) and neglecting the other ones. In the 



resulting equations only the zero-momentum as well as 
the k and — k atomic and molecular modes are present 
which allows to solve the coupled linear operator equa- 
tions numerically. Their asymptotics (t — > oo) may be 
investigated analytically though, yielding the following 
expectations values for the number operators: 
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where A = y / 2Nil 2 — £ 2 n 4 . A is real for all modes with 
n < ^/2NQ 2 /£ 2 , which sets the number of modes whose 
population grows in time. As typically both f2 and £ 
are of the order of 10 - 10 2 Hz while N ~ 10 5 - 10 6 
f7|, ^0|JT5| , corrections to the 2-mode model are divergent 
for many low-lying states. With the parameters used it 
is only for the 11-th and higher excited atomic modes 
that the quantum corrections are small and oscillatory 
(imaginary A). Therefore, one is not allowed to exclude 
excited modes (k ^ 0) from a theoretical model Q . 
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FIG. 1. Populations of the atomic (solid line) and molec- 
ular (dotted line) condensates calculated within a 2-mode 
model (see Eqs. (H)) and of the quantum corrections (dashed 
line) in the 1st excited atomic mode (see Eqs. (^)). The pa- 
rameters are ft =18.432 Hz, £=71.373 Hz and N = 10 5 . 

Fig. |l| presents an example of the run-away correction to 
the 2-mode model obtained numerically together with the 
classical source terms cvo(i) and Po(t) (note: in all plots 
the total populations are normalized by the total number 
of particles N). From Fig. [i] one immediately concludes 
that the 2-mode solution (|5|) is physically invalid for times 
larger than ~ 0.6 ms. To cure this problem, from now on, 
we will use a multi-mode model in the form of Eqs. (M). 

Solution of the nonlinear operator equations (3) 
presents an extremely difficult task. A semiclassical ap- 
proximation, however, is well justified for all except ex- 
tremely low temperatures. Therefore, we replace all oper- 
ators by c-number complex amplitudes. From the view- 
point of the Bogolubov method [12], such an approach 
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is legitimate as indeed many modes are macroscopically 
populated (i.e. their occupation is greater than quantum 
fluctuations) . This way we are left with a set of nonlinear 
differential equations which must be solved numerically. 

The first observation in the multi-mode model is that 
if one starts from a pure atomic condensate (the k = 
mode), the 2-mode dynamics is recovered. However, even 
a very small occupation of excited atomic or molecular 
modes results in the dynamics beyond the 2-mode ap- 
proach. Such a behavior resembles an initiation of super- 
fluorescence where quantum, not thermal, fluctuations 
play a role. In a typical experiment, roughly 85% of the 
total number of atoms populate the atomic condensate 
whereas the rest of them is thermally distributed over 
excited atomic modes and we mimic such a situation in 
the initial conditions of our model jl7) . All the molecu- 
lar modes are initially unpopulated. Each atomic mode 
is assigned an initial, randomly chosen, phase. Any sub- 
sequent dynamics depends on the initial phases and, in a 
sense, a single simulation describes a single experimental 
realization. Values of the parameters in the model are 
fl =18.432 Hz, £ =71.373 Hz, N = 10 5 and g =0.018 
Hz (the atomic mass and the scattering length are those 
of 87 Rb and the size of the box is equal to the Thomas- 
Fermi radius of a condensate of N atoms in a trap with 
frequency of lu — 2ir 80 Hz - see pp| , p!5| ). Sample re- 
sults for the model with 2622 modes (maximum atomic 
and molecular excitation rii = —5, . . . , 5, i = x, y, z) are 
presented in Fig. |2[ In fact the results stabilize for the 
number of modes exceeding 1000. 
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FIG. 2. Populations of atomic (thick solid line) and molec- 
ular (thick dashed line) condensates (the k = modes) and 
the sums of populations of all excited (k ^ 0) atomic (thin 
solid line) and molecular (thin dashed line) modes in the 
model with 2622 states. 

A striking discovery is that after a relatively short time 
both the atomic and molecular condensates (the k = 
modes) are completely depleted and all particles occupy 
excited modes in roughly equal proportions. For the used 
(typical) parameters, only one oscillation in the conden- 



sate population survives - see Fig. |3|. 
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FIG. 3. 2-mode (dashed line) vs. multi-mode (solid line) 
dynamics of an atomic condensate population. 

This clearly indicates that if one wants to convert an 
atomic condensate to its molecular counterpart, it is nec- 
essary to precisely tailor the length of the coupling pulse 
so that the maximum of the molecular k = mode is 
picked. If one investigates the condensate (k = 0) pop- 
ulation and a sum of populations of all excited modes 
(k 7^ 0) solely (like in Fig. ||), their dependence on the 
initial phases is negligible in a sufficiently big model (we 
checked this fact by setting different initial phases in our 
calculations). The effect is due to self-averaging caused 
by very many random phases present in the system. How- 
ever, it can be observed in the time evolution of single 
modes. We emphasize that, contrary to the previous 
treatments j^-|ll],|l5| , the effective losses from the con- 
densates are completely due to a Hamiltonian evolution 
and not because of any phenomenologically introduced 
loss processes. In other words, due to an external weak 
coupling of many degrees of freedom, the system is ef- 
fectively heated and its final equilibrium state certainly 
does not result from any T — dynamics. The inclusion 
of atom-molecule and molecule-molecule collisions would 
randomize the still coherent dynamics leading to an anal- 
ogous and even more pronounced effect. Remarkably, the 
character of the dynamics does not depend on the total 
number of particles - our simulations for ./V = 2 ■ 10 3 and 
TV = 5 • 10 5 still show only one oscillation in the atomic 
condensate population on a slightly altered (longer for 
smaller N) time scale. The amplitude of this oscillation 
is bigger for smaller N indicating a small increase in the 
molecular condensate production (~ 30% for N = 2 • 10 3 
and ~ 23% for N = 5 • 10 5 instead of - 25% for N = 10 5 
as seen in Fig. |J). 

In order to recover the standard oscillatory dynamics 
of the 2-mode limit P-pH within a multi-mode model, 
one needs to detune the excited modes. The latter can 
be achieved by decreasing the size of the box (and so 
enlarging a spacing between the excited modes) while 
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keeping the total density fixed (in order not to alter the 
scattering and coupling parameters). The results for the 
box whose volume is 15 3 times smaller are presented in 
Fig. ^ (the appropriately decreased number of particles 
is N = 10 5 /15 3 ^ 30). 
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FIG. 4. Populations of atomic (thick solid line) and molec- 
ular (thick dotted line) condensates (the k = modes) and 
of the sums of all excited (k 7^ 0) atomic (thin solid line) and 
molecular (thin dotted line) modes in very small box (see text 
for details) - recovery of the 2-mode dynamics. 



To summarize, we point out the invalidity of a 2-mode 
model for a description of present experiments with cou- 
pled atomic-molecular systems in the Bose-Einstein con- 
densation regime. We find it necessary to employ a multi- 
mode approach and within it we observe a complete de- 
pletion of both the atomic and molecular condensates on 
a short time scale. For typical parameters, only one giant 
oscillation between the two condensates is present. Thus, 
a new destructive mechanism in the system is pointed out 
which should be taken into account while planning future 
experiments. An effectively 2-mode dynamics is recov- 
ered under special conditions of large detunings. The 
system presents an interesting example of a quantum co- 
herent dynamical evolution effectively randomized by the 
coupling of many degrees of freedom. 

In this Letter we presented the calculations for the box 
rather than the harmonic oscillator potential. We have 
chosen the box size such that the first excitation energy 
is of the order of the level spacing in the experiment of 
Heinzen [Q] . Hence, we may expect that for the harmonic 
oscillator the time of coherent evolution would be even 
shorter, because of the quadratic vs. linear dependence 
of excitation energies in the respective potentials. 
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At the time of writing this manuscript we found out 
that the problem of going beyond the 2-mode dynam- 
ics had been also addressed in M. Holland, J. Park, and 



R. Walser, cond-mat /0005062 The authors point out an 
effect of condensate depletion due to a thermal cloud. 
However, in their treatment they neglect all collisions. 
We do recover the 2-mode dynamics in the semiclassical 
approximation starting from the pure atomic condensate 
(very low temperature). However, as we see in Fig. [j], even 
in the case of a 100% initial occupation of the atomic 
condensate, the quantum initiation makes the 2-mode 
model invalid after ~ 0.6 ms. 
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